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Abstract — Typical adaptive controllers are restricted to using 
a specific update law to generate parameter estimates. This 
paper investigates the possibility of using any exponential 
parameter estimator with an adaptive controller such that the 
system tracks a desired trajectory. The goal is to provide 
flexibility in choosing any update law suitable for a given 
application. The development relies on a previously developed 
concept of controller/update law modularity in the adaptive 
control literature, and the use of a converse Lyapunov-like 
theorem. Stability analysis is presented to derive gain conditions 
under which this is possible, and inferences are made about the 
tracking error performance. The development is based on a 
class of Euler-Lagrange systems that are used to model various 
engineering systems including space robots and manipulators. 

I. Introduction 

Control laws with adaptive feedforward terms have been 
widely used to compensate for the linear-in-the-parameter 
(LP) uncertainty in nonlinear systems. The update law or 
a parameter estimator is an integral part of an adaptive 
controller that generates parameter estimates used by the 
adaptive controller. The convergence properties of the es- 
timator directly impact the tracking error performance, even 
though asymptotic tracking may still be possible without 
the estimates converging to their true values. Lyapunov- 
based techniques are commonly used to design and analyze 
an adaptive control law and its corresponding parameter 
update law, where the controller and the update law are 
designed in conjunction with one another such that the update 
law cancels certain cross terms in the stability analysis to 
make the overall time-derivative of the Lyapunov function 
negative. This restrains the design of the update law to 
one specific form, which typically is a gradient update law. 
Gradient update laws are notorious for exhibiting slower 
parameter convergence [1] — [3], which could lead to a slower 
transient performance of the tracking error in comparison 
to other possible adaptive update laws (e.g., least-squares 
or recursive least-squares update law [1], [3]— [5]). In order 
to overcome this issue of slow convergence, several results 
have been developed in the literature that aim to modify or 
augment the typical tracking error-based gradient update law. 
For example, composite adaptive update laws [1], [2], [6]; 
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prediction error (a predicted value of the difference between 
the actual parameter and its estimate)-based update laws [7] — 
[11]; and various least-squares update laws [3], [5], [12], In 
addition, several results have discussed the robustness and 
exponential convergence properties of estimation schemes 
[1], [4], [5], [13], [14], However, the use of these update laws 
is still contingent upon their ability to cancel cross terms in 
the Lyapunov-based stability analysis. That is, only a specific 
type of update law can be used with the adaptive controller. 
This raises the question if we can have flexibility in choosing 
a parameter update law best suited for a given application to 
potentially achieve a better performance. 

Intuitively, if a parameter estimator has an exponential 
convergence, it could be combined with the adaptive con- 
troller to achieve tracking. However, as mentioned earlier, 
stability cannot be guaranteed unless it is designed to cancel 
certain cross terms. The result in this paper is motivated by 
the desire to investigate if a formal proof can be presented for 
guaranteeing stability and tracking performance if any such 
exponential parameter estimator is used with the adaptive 
controller. In the existing literature, there are various results 
on a class of modular adaptive controllers that provide 
such controller/update law separation fcf. [7], [9]— [1 1]). One 
approach is to use nonlinear damping [8], [15] to yield 
a strong controller that is input-to-state stable (ISS) with 
respect to the parameter estimation error. This enables the use 
of any update law that yields bounded parameter estimates. 
A brief survey of nonlinear damping-based modular adaptive 
control results is provided in [7], Another approach was 
investigated in [16] by designing a robust control law using 
the RISE (robust integral of the sign of the error [17]) 
feedback method, such that a generic form of an update law 
in terms of tracking errors could be used. 

In this paper, modularity in the update law (i.e., con- 
troller/identifier separation), and tracking is achieved in 
three steps. First, for a system which is linear in uncertain 
constant parameters, an adaptive controller is developed and 
augmented with a nonlinear damping term. This makes the 
system ISS with respect to the estimation error. That is, for 
a bounded parameter estimation error, the tracking error is 
bounded. Second, a converse Lyapunov-like theorem is pre- 
sented, and existence of a Lyapunov-like function is shown 
for systems with trajectories that are exponentially decaying 
to a small bound. Third, any parameter estimator that gives a 
parameter estimation error exponentially decaying to a small 
bound is used with the adaptive controller. Based on the 
second step, the existence of a Lyapunov-like function for 
such a parameter estimator is guaranteed. This combined 



with the first step is used to prove a tracking result under 
some sufficient gain conditions such that, if the parameter 
estimation error is uniformly ultimately bounded (UUB), 
then the tracking error is also UUB. 

Typical adaptive control designs yield asymptotic tracking 
without any tracking performance guarantees. In contrast, 
the current result proves only UUB tracking, but guarantees 
performance bounds on the tracking error while allowing the 
use of a broad class of exponential estimators. Moreover, if 
the parameter estimator has an exponential convergence to 
zero, the tracking errors also converge to zero exponentially. 
Compared to the previous results on controller/update law 
modularity, the current result can be considered as an alterna- 
tive choice by the control design engineers. The result in this 
paper achieves modularity in the controller/update law using 
ISS arguments similar to [7]— [1 1] by including a nonlinear 
damping term in the torque control input. However, in [7]- 
[11], in order to prove tracking, additional assumptions about 
the update law itself had to be made. For example, it needed 
to be driven by a prediction error, which must be square 
integrable; the update law must be bounded, etc. This in 
turn required extra terms in the control law. Moreover, the 
estimator must also ensure that the estimate of the inertia 
matrix be positive definite. In most cases, this required the 
use of projection algorithms, thus leading to an increased 
complexity, and even though asymptotic tracking was proved, 
no bounds on the tracking error transient could be derived. 
In contrast, the result in this paper only assumes that the 
parameter identification scheme gives an estimation error that 
exponentially decreases to a small bound, and no additional 
assumptions are required. 

Since the objective is to design an adaptive control law 
that utilizes parameter estimates generated by any estimator 
with exponential convergence, design of a specific estimator 
is not considered in this work. It should be noted that the 
construction of any such practical estimator may require 
additional inputs, such as bounds on measurement noise, 
bound on estimation error, etc. [1], [4], [5], [13], [14]. 

II. Dynamic Model and Properties 

The class of nonlinear dynamic systems considered in 
this paper can be described by the following Euler-Lagrange 
formulation, which can be used to model a large class of 
engineering systems such as space robots and manipulators 
[18], [19]: 

M(q)q + V m (q,q)q + G(q) + F (q) = r(f). (1) 

In (1), M(q) £ R" xn denotes the inertia matrix, V m (q, q) £ 
R" x " denotes the centripetal-Coriolis matrix, G(q) £ R™ 
denotes the gravity vector, F (q) £ R n denotes friction, 
r(f) £ R™ represents the torque input control vector, and 
q(t), q(t), q(t) £ R" denote the generalized position, veloc- 
ity, and acceleration vectors, respectively. The subsequent 
development is based on the assumption that q(t ) and q(t.) 
are measurable and that M{q), V m {q , q), G(q), and F (q) are 
unknown. Throughout the paper | • | denotes the absolute value 
of the scalar argument, ||-|| denotes the standard Euclidean 


norm for a vector or the induced norm for a matrix. The 
following properties will be exploited in the subsequent 
development. 

Property 1: The inertia matrix M(q) is symmetric, posi- 
tive definite. Therefore, it satisfies the following inequality 

Vg(*U(t) GR” [18], [19]: 

mi ||£|| 2 < i r M{q)( s < m 2 ||£|| 2 ( 2 ) 

where mi, m 2 G R are known positive constants. 

Property 2: The time derivative of inertia matrix M (q) and 
the centripetal-Corilois matrix V m (q, q) satisfy the following 
skew-symmetry property Vq(t),q(t),£(t) £ R” [18], [19]: 

£ T (^M(q)-V m (q,q)^J £ = 0. (3) 

III. Control Objective and Development 

The objective is to design an adaptive controller with 
the flexibility of using any parameter estimator that has an 
exponential parameter convergence to a small error bound, 
such that the system tracks a desired time-varying trajectory 
q,i{t) despite uncertainties in the dynamic model. To quantify 
this objective, a position tracking error, denoted by ei (t) £ 
R", is defined as 

ei = q,i - q ■ (4) 

In addition, to facilitate the subsequent analysis, a filtered 
tracking error [18], 62 (f) G R n , is defined as 

e 2 = ei + Qiei (5) 

where a\ £ R denotes a positive constant. Differentiating 
(5) with respect to time, and premultiplying the resulting 
expression with M (q) yields 

M (q) e 2 = M (q) q d - M (q) q + M (q) a.\e\. (6) 

Using (1), (4), (5) in ( 6 ), we get 

M (q)e 2 = Y6- r(t) - V m (q, q)e 2 (7) 

where the following property is used: 

Property 3: A part of the dynamics in (1) can be linearly 
parameterized as [18], [ 20 ] 

Y8 = M(q)q d + V m (q,q)(q d + a 1 e 1 ) +G(q) 

+F (q) + M (q) a 1 e 1 (8) 

where 0 G R 7 ' contains the constant unknown system param- 
eters (e.g., component masses, moments of inertia, damping 
and friction coefficients, etc.), and Y (q, q , q d , q d ) £ R" xp is 
a regression matrix that contains known nonlinear functions 
of the known generalized coordinates q(t), q(t) and desired 
trajectories. 

Based on the open-loop error system in (7), the torque 
control input is designed as 

r(f) = Y6 (f) + fcie 2 (f) + k n \\Y || 2 e 2 (f) (9) 

where ki , k n £ R are positive constants, 6 (f ) G R p is 
an adaptive estimate of 8. Nonlinear damping terms, similar 
to k n || y || 2 e 2 (f) in the control input, have been previously 



used in results such as [8], [10]. Substituting (9) in to (7), 
we get the following closed-loop error system: 

M (q) e 2 = Yd-kie 2 - k n ||U|| 2 e 2 - V m (q,q)e 2 (10) 

where 9 (t) £ R p is the parameter estimation error defined 
as 

d{t) = e-e(t). 

IV. INPUT-TO-STATE STABILITY 

In this section, we show the ISS property of the system 
with respect to the parameter estimation error 9 (t). This 
allows controller/update law modularity, and thus, any update 
law that generates bounded parameter estimates can be used 
with the adaptive controller. We now have the theorem: 
Theorem 1: The closed-loop error system in (10) is input- 
to-state stable with respect to 9 ( t ), i.e., if 9 ( t ) £ Coo, 
then ||ei (t)\\ € Coo- 

Proof: Define a Lyapunov function candidate 

v iss = ^2 M {q) e 2 . (11) 

Using (2), Vi SS (t) can be lower and upper bounded as 



Using the fact that a 2 < b 2 + c 2 => a < b + c for all positive 
scalars a, b , c, (17) can be expressed as 


— ||e 2 (0)|| e (-A ” 5t) 

mi 


= sup 9 (t) 

2 VTOi A i s s k n 0<r<£ 


The bound in (18) indicates input-to-state stability (ISS) 
(cf. Def 4.7, pp. 175 [21]). That is, if sup 0<T<i 9 (r) is 
bounded, then ||e 2 (t ) || is bounded, and therefore, from (5), 
the tracking error ||ei (f)|| is bounded. □ 

Any update law that guarantees a bounded 9 can be 

used. Further, if 9 ( t ) is UUB, then subsequently it will 
be shown (Theorem 3) that the tracking error is UUB. First 
we prove the following converse Lyapunov-like theorem. 


\m 1 ||e 2 || 2 < Viss < ^m 2 ||e 2 || 2 . (12) 

Taking the time derivative of (1 1), and substituting from (10), 
we get 

Viss = e% [Y9- he 2 - k n ||U|| 2 e 2 - V m (q, g)e 2 ](13) 
+ \e^M (q) e 2 . 

Using the skew-symmetry property (3), the expression in (13) 
can be upper bounded as 

Viss < ~ki ||e 2 || 2 + ||e 2 || ||y|| ||#|| — k n ||e 2 || 2 ||V|| 2 . 

Completing the squares for the last two terms 

Vu. < -h ||e 2 || 2 -k n (j|F|| || e 2 1| - +^~ ||0|| 2 

therefore, 

V S s < — k\ ||e 2 || + — — 9 . (14) 

4/C n 

Using (12) in (14), the upper bound on V)ss (t) can be 
expressed as 

Viss ^ — 2A issViss + 77“ 9 (15) 

4/C n 

where 

-h. 

/ 'l SS — 

m 2 

After solving the differential inequality in (15), we get 

Viss (t) < Viss (o) eV 2 ^^ + 1 sup 9 ( T ) . 

OAi ss rC n 0<r<t 

(16) 


V. Converse Lyapunov-like Theorem 

The following can be proved directly by extending the 
result in Theorem 4.14 of [21], 

Theorem 2: Consider a nonlinear system 

x(t) = f(t,x) (19) 

where / : [0, oo) x D -A R" is continuously differentiable, 
D = {x £ R n | ||a;|| < r} for some r, and the Jacobian 
matrix [df / dx\ is uniformly bounded on D. Let x = 0 be 
an equilibrium point for the system in (19). Let k, A, e, 
and ro be positive constants with ro < r/k. Let Dq = {x £ 
R” | ||a:|| < r 0 }. Assume that V x (to) £ D 0 , V t > to > 0, 
the trajectories of the system in (19) satisfy 

\\x (f)|| < k ||x (t 0 ) || e _A(t_to) + e. 

Then, there exists a function V : [0, oo) x Do —> R that 
satisfies the inequalities 

ci ||a;|| 2 < V (t,x) < c 2 ||a:|| 2 + c 3 ||x|| + c 4 

dV dV 9 

~dt + ~dx ^ X ^ + ° 6 11*11 + ° 7 

for some positive constants ci, c 2 , . . . , C7. Moreover, if r = 
00 and the origin is globally uniformly ultimately bounded, 
then V (t. x) is defined and satisfies the aforementioned in- 
equalities on R". Furthermore, if the system is autonomous, 
V can be chosen independent of t. 

Proof: Let (f> (r; t, x) denote the solution of the system start- 
ing at (t, x); i.e., <j) (t; t, x) = x. For all x £ D 0 , </> (r; t, x) £ 
D for all t >t. Let 

pt -\-8 

V (t, x)= 4 > t (t; t, x) (j> (r; t, x) dr (20) 



where S is a positive constant to be chosen. Due to expo- 
nentially decaying bound on the trajectories, we have 


rt+8 

v (t, x) = J \\<j> (r; t, x)|| 2 dr 


< J {k 2 \\<t>(t-,t,x)\\ 2 e 2A(t t] 

+2 ek || cj>(t\ t, a?) || e _A ( r_t ) + e 2 )dr 
Using the fact (f> (f; t, x) = x, 

V 0 1 , x) < ^ (1 - e~ 2XS ) \\x\\ 2 + 2 -^ (1 - e"^) |M|+ £ 2 A 


we get 

^ + ^/(f,a;) = 4> T (t + 8;t,x)<j>(t + 5;t,x) - ||a:|| 2 

< (k ||x|| e~ xs + e) 2 — ||x|| 2 

< -(l-k 2 e~ 2XS )\\x\\ 2 
+2 eke~ xs ||a;|| + e 2 . 

Therefore, the second inequality of the theorem is satisfied 
with 

c 5 = 1 — k 2 e~ 2XS , c 6 = 2 eke~ xs , c 7 = e 2 . (27) 


Since, the Jacobian matrix [df/dx] is uniformly bounded on 
D, let 

df 

— (t,x) < L, V x € D. 


\\f(t,x)\\<L\\x\\ (22) 

and <fi (r; t, x) satisfies the lower bound (see Appendix, 
Lemma 1) 

110 ( r ! t, at) || 2 > ||x|| 2 e _2i(r_t) . (23) 

Hence, using (23) in (20), V (£, x) satisfies the lower bound 

rt+8 1 

V(t,x)> j e- 2L ^dT\\x\\ 2 = -(l-e- 2LS )\\x\\ 2 . 

(24) 

Thus, from (21) and (24), V (f, x) satisfies the first inequality 
of the theorem with 

C = ft = Yx h - e "“‘) (25) 

c 3 = c i =E 2 S. 

To calculate the derivative of V (f , x) along the system 
trajectories, the following sensitivity functions are defined 

d d 

0t 0; t, x) = — 0 (r; t, x) ; (j) x (r; t, x) = — 0 (r; t, x) . 


dV dV 

— + — /(f,a;) = c/) T (t, + 6;t,x)(l)(t + S;t,x) 

—<j) T (f; t, x) </> (f; t, x) 

rt+8 

+ 2cj) T (r;t,x) 4> t (r;t,x) dr 
rt+5 

+ 2+ (r; t, x) 4> x (r; t, x) dr 

xf(t,x) 

= (f) T (t + (5; t, x) 0 (t + J; t, x) 

rt+5 

+ J 2+ (r; t, x) [<f> t (t; t , x) 

+0* (r; t, x) f ( t , x)]dr - ||a;|| 2 . 
Using the fact that (see Appendix, Lemma 2) 

0t (t; t, x) -(- 0 X (r; t, x) f (t, x) = 0, Vr>f (26) 


VI. Parameter Estimation and Tracking 

Theorem 3: Suppose there exists a parameter estimator such 
that the parameter estimation error 0 ( t ) satisfies 

0(f) < fc 0 (0) e~ xt + e (28) 

then, the control law (9) in conjunction with the parame- 
ter estimator yields globally uniformly ultimately bounded 
(GUUB) tracking stability for the system in (1) and all 
closed-loop signals are bounded. 

Proof: Theorem 2 implies that there exists a function 
V c if (f:, 0) such that 

,, II 2 || ,| 2 || ,| 

ci i < V c if < C 2 0 + c 3 0 + C 4 (29) 


Define 


Vdf < -C 5 0 +c 6 0 +c 7 . 


V a — Vi ss + Vdf 


where Vi ss was defined in (11), which satisfies the following 
lower and upper bounds: 

V s > ||e 2 || 2 + ci 0 (32) 

V s < ^to 2 ||e 2 || 2 +c 2 0 +C3 0 +C4. (33) 

Taking the time derivative of (31), and using the results of 
Theorems 1 and 2 

V s <-ki\\e 2 \\ + —— 0 — C5 0 +C6 0 +c 7 . (34) 

4/C n 

Using the upper bound of (33) in (34) 


Vs ^ 2A iss^s ( C 5 e Z\i ss C2 


H - (2A2SS C3 ~b Cq) j j ^|| ~b 2Xi ss C4 + C7. ( 35 ) 

Let 'yi — C5 2AjssC2 ’ 3^2 — 2A^ SiS C 3 ~b C6, ^3 — 
2 A^ s<s C 4 + C7, where 71, 72, 73 are positive constants such 
that the following sufficient conditions should be satisfied 

kn> +c 5 -2\ issC2 y Xiss< 2+y 


(36) 



Completing the squares for the second and the third terms 
in (35) 

V s < -2X iss V s + 7 (37) 


where the positive constant 7 is defined as 

7? 

7 = ^~ +73- (38) 

471 

After solving the differential inequality of (37), we get 


V a (t) < V. (0) e 

From (32), we have 


-2A, 


7 


2A, 


(1- 


-2A i 


(39) 


^1 l|e 2 || 2 < ^mr ||e 2 || 2 + Ci 


< 


Then, the following bound on the filtered tracking error e 2 (t) 
can be obtained: 


|„ /J.M|2 2C S (0) —2X iss t 

e 2 (t) < e 

TOi 


7 


miAi 


Therefore, 


|e 2 (f)||< 


I 2V S (0) x _ t 

mi 


mi A,;, 


(40) 


(41) 


Using (5), a similar bound for ||ei (t)|| can be derived, and 
hence, q(t) , q (t) are bounded and the control torque (9) is 
bounded. □ 

Remark 1: Substituting for c 2 and C5 from (25) and (27) in 
to the gain condition for Xi SS in (36), we have 


A iss Y A 


(l — k 2 e 2A<5 ) 
k 2 (1 — e~ 2XS ) ' 


(42) 


The expression on the RHS is of the order of O (A). Thus, the 
condition in (42) indicates that the decay rate of the tracking 
error is less than that of the estimator. In other words, the 
parameter estimator should have a fast convergence in order 
to ensure a good tracking performance. Similarly, from (40) 
using (38), we have 


7 


1 


miAi, 


miXi, 

1 

miXi, 


'± 

,47i 


■73 


(2Xi ss C3 + Cq) 2 


+ 


1 


miAi, 


^4(c 5 2Ai SS c 2 1/4 k n ) 

(2A i ss c 4 + c 7 ) . 


Substituting Ci, c 2 , . . . , c 7 from (25) and (27), 


VII. Simulation Study 


In order to validate the theoretical development, a com- 
puter simulation was performed on a space robot system. 
Control of space robots is important for various space appli- 
cations such as material transport, cooperative autonomous 
construction, maintenance of space stations and satellites, 
etc. Specifically, we consider a space robot consisting of a 
base (mass mo), and a two-DOF revolute manipulator with 
equal link lengths l and point masses mi and m 2 at the end 
of each link. Although the proposed method can control the 
full system, for simplicity the proposed adaptive controller is 
applied only for the motion control of the arm. It is assumed 
that the attitude of the base can be effectively controlled by 
other control methods. Neglecting the gravity term in (1), the 
system dynamics can be derived in a manner similar to [22] 
so that M = P 1 R 1 +P 2 R 2 +P 3 R 3 , where pi = m ^ 1 , p 2 = 
"ra™ 2 , P 3 = 77, where m c = mo + mi + m 2 , and 


Z 2 O' 


0 0 


i?2 


z 2 r- 

1 2 r- 


where c 2 = cos (q 2 ). V m = p 3 R 

— 2 l 2 s 2 q -2 ~l 2 s 2 q 2 

l 2 s 2 qi 0 


i?4 = 


where s 2 = sin (g 2 ) ■ 


A continuous friction model similar to [23] is used, 
F = 7ii? 5 + 7 2 i? 6 + 7 3 <j, where R 5 = tanh(O.lq) - 
tanh(0.2g), R 6 = tanh(0.15g). It is assumed that the 
link length l is known, while the masses mo, mi, m 2 
(equivalently, pi, p 2 , p 3 ) and the friction coefficients 
71, 72, 73 are unknown. Defining = q d + aaei, C2 = 
id + aqei, and comparing with (8), we have Y = 
[R 1 C 1 ^2(2 (-R3C1 + -R4C2) R 5 Rq q] and the un- 

known parameter 9 = [pi p 2 p 3 71 7 2 73] . An 
arbitrary parameter estimator is constructed such that the 
bound in (28) is satisfied. The parameter estimates are 
initialized to half of their true values and the resulting 
parameter estimation error is shown in Fig. 1. The adaptive 
control law in (9) with the control gains chosen by trial and 
error is employed to track a desired sinusoidal trajectory. The 
position tracking error converges to a small bound (shown 
in Fig. 2) as proved by the stability analysis. 


Fig. 1 . 9 (t) exponentially converges to a small bound. 


7 

m 4 A i s 


4 e 2 k 2 (2X iss (l - e xs ) + Ae x5 ) 2 
(A - Xk 2 e~ 2XS - X iss k 2 (1 - e~ 2XS ) - 

1 ^2 


x 




which is of the order of O (e 2 ), and therefore, from (41), 
the ultimate bound on ||e 2 (t)\\ is O (e). 

Remark 2: If the estimation error 9 (t) in (28) exponentially 
converges to zero (i.e., e = 0), then from (43), the tracking 
error also exponentially converges to a zero steady-state. 



Fig. 2. Position tracking error. 



VIII. Conclusion and Future Work 

Using a converse Lyapunov-like theorem, sufficient gain 
conditions were derived under which any exponential pa- 
rameter estimator can be used with an adaptive controller to 
achieve tracking with a performance bound on the tracking 
error. The adaptive controller is augmented with a nonlinear 
damping term which makes the closed-loop system input-to- 
state stable enabling controller/update law separation. The 
current development is based on Euler-Lagrange dynamic 
systems. Future work will focus on extending the result to 
more general systems with parametric uncertainty. 
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IX. Appendix 

Lemma 1: In (23), \\(j) (r; t, x)\\ 2 > ||x|| 2 e _2i ( T_t ). 

Proof: From (19) and (22) 


j f [x T (t) x (t)] 


2x t (t) x ( t ) | 

2x T (t) f (t,x)\ < 2L\\x\\ 2 . 


Therefore, 

-2L ||x|| 2 < ~ [x T (: t ) x (t)] < 2 L ||x|| 2 . (44) 


Solving the differential inequality in (44) such that 4> (r; t, x) 
denotes the solution of the system starting at (t, x) (i.e., 

<f> (f; t, x) = x) 


\\x\\e~ L ^ <U{r-,t,x)\\< |N|e L ^. 


Hence, (23) holds. □ 

Lemma 2: In (26), x) + (f> x (T;t,x) f (t,x) = 

0, V r > t. 

Proof: The solution <f>(r m ,t,x) starting at (t, x) can be 
expressed as 

0(r;f,x) = <j> (f; t,x) + J f (s,(j>(s))ds 
= X + J t f {sA(s))ds. 

Therefore, 

fix 

(j)t (r; t, x) = — - f (f, <t> (t)) = -/ (t, x) (45) 

and „ 

OX 

= — = 1. (46) 

From (45) and (46), it can be verified that (26) holds. □ 


